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Parallel  Computation  of  Three-Dimensional, 
Unsteady  Wake  Flows  Using  Vortex  Methods 


Technical  Objectives 

The  objectives  of  the  research  were  to  develop  computational  techniques  for  parallel  computers  capable 
of  simulating  complex,  three-dimensional,  unsteady  wake  flows  at  high  Reynolds  numbers  and  to  apply 
these  techniques  to  such  flows  as  those  about  maneuvering  three-dimensional  bodies.  A  Lagranpan  vortex 
method  was  used  as  a  basis  for  the  numerical  scheme.  In  vortex  methods  the  only  grid  required  is  that 
to  mesh  the  surface  of  the  body.  In  addition,  the  computational  elements  containing  vorticity  move  with 
the  local  fluid  velocity  and  are  required  only  in  the  boundary  layers  and  wake  of  a  bluff  body  flow.  This 
technique  is  therefore  a  natural  candidate  for  the  simulation  of  wake  flows  in  that  the  fluid  containing 
vorticity  for  these  flows  occupies  only  a  small  fraction  of  the  total  fluid  volume  of  interest. 

Technical  Approach 

A  robust  vortex  particle  method  for  three-dimensional  flows  was  developed,  including  a  new  treatment 
of  viscous  effects  to  allow  for  the  accurate  treatment  of  boundary  layer  mechanics  including  unsteady 
separation.  In  addition,  to  allow  for  ambitious,  three-dimensional  applications,  fast  algorithms  were 
developed  for  both  the  vortex  particle  method  and  the  panel  method  (used  to  satisfy  the  no-slip  boundary 
condition  on  the  body  surface).  These  algorithms  were  implemented  on  the  Intel  Delta  supercomputer  at 
Caltech.  This  576  mode  machine  has  528  Intel  i860  numerical  processors  giving  a  total  computing  capability 
of  over  30  Gflops  peak  performance  with  8  Gbytes  of  memory. 

Accomplishments 

We  carefully  considered  three  possible  choices  for  a  fast  algorithm  with  regard  to  their  speed,  accuracy 
and  implementation  on  a  parallel  computer.  The  various  required  translation  operators  were  developed 
along  with  the  necessary  advanced  error  estimation  techniques. 

By  working  in  close  collaboration  with  J.  K.  Salmon  (Research  Fellow,  Caltech  Physics  Department) 
we  modified  several  versions  of  fast  parallel  tree  codes  developed  for  large  scale  3D  gravitational  N-body 
problems  into  codes  for  3D  fluid  dynamical  N-vortex  particle  problems.  On  the  Intel  Delta  with  512 
processors,  inviscid  and  viscous  computations  involving  N  =  10^  to  10^  3D  vortex  particles  were  earned 
out  successfully.  Timings  and  scaling  analyses  were  completed,  covenng  the  range:  P  -  2"  proce^ore 
(n  =  0, 1 , 2,  ....9)  and  W  =  10^  -  10^  vortex  particles.  They  show  that  the  CPU  time,  T ,  scales  as  T  a 
and  that  the  parallel  overhead  (obtained  by  measuring  the  difference  between  the  P  =  512  curve  and  the 
extrapolation  of  the  P  =  1  curve  for  large  N)  is  in  the  neighborhood  of  20%. 

In  addition,  the  fast  parallel  tree  code  for  vortex  particles  was  itself  modified  into  a  code  for  boundary 
element  methods  (BEM)  (also  called  panel  methods).  In  its  present  implementation,  panels  that  are  infinitely 
thin  vortex  sheets  are  used.  In  the  context  of  the  present  application,  the  vorticity  represented  by  these 
panels  is  that  required  to  satisfy  the  no-slip  boundary  condition  at  the  body  surface  by  creating  the  proper 
amount  of  boundary  vorticity  (i.e.,  "bound"  vorticity)  at  every  time  step.  This  fast  parallel  tree  code  for 
M-panel  problems  has  been  used  to  compute  potential  flows  past  simple  geometries  such  as  the  sphere  an 
the  ellipsoid  at  angle  of  attack.  Timings  were  been  carried  out,  spanning  again  the  ranges  ^  =  2  processors 
(n  =  0, 1,  ....,9)  and  M  =  10^  -  10®  vortex  sheet  panels,  with  the  same  type  of  scaling  as  obtained  for  the 
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The  Ph.D.  thesis  work  of  P.  Koumoutsakos  which  was  completed  3  September  1992  gave  us  an  excellent 
understanding  of  the  creation  and  unsteady  transport  and  diffusion  of  vorticity  near  a  solid  boundary  vis-a-vis 
2D,  viscous  vortex  methods. 

Significance  of  Accomplishments 

The  development  and  implementation  of  the  fast  parallel  tree  code  for  3D  vortex  particles  on  the 
Intel  Delta  constitutes  a  major  building  block  of  the  present  program,  and,  along  with  computations  with 
AT  =  10^  -  10^  3D  vortex  particles,  represents  a  major  accomplishment  in  the  field  of  vortex  dynamics. 
In  a  demonstration  computation  we  studied  the  time-evolution  of  an  initially  spherical  vortex  sheet  with 
N  =10^.  This  proved  very  interesting  and  detailed:  3D  instabilities,  vortex  structure  interactions,  etc. 

The  development  and  implementation  of  the  fast  parallel  tree  code  for  3D  boundary  element  methods 
on  the  Intel  Delta  also  constitutes  a  major  building  block  of  the  program,  and  a  significant  accomplishment. 
Whereas  the  state  of  the  art  in  potential  flow  computations  is  usually  in  the  range  Af  =  10^  -  10'^  panels,  we 
can  now  contemplate  computing  potential  flows  with  M  =  10^  —  10^  panels.  This  fast  panel  code  should  be 
of  great  interest,  not  only  in  fluid  mechanics:  potential  flow  past  bodies  of  very  high  complexity,  free  surface 
flows  (e.g.,  marine  engineering,  wave  dynamics,  ship  +  wave  dynamics),  bubble  flows  (e.g.,  chemical 
engineering),  etc.;  but  also  possibly  in  other  disciplines:  electrostatics,  radar,  antennas,  magnetostatics, 
acoustics,  in  fact  most  boundary  integral  problems  expressed  in  terms  of  a  Fredholm  integral  equation  of 
the  second  kind  and  solved  using  a  boundary  element  (=  panel)  method. 
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Anthony  Leonard,  professor  of  aeronautics,  is  the  PI.  Petros  Koumoutsakos,  a  postdoctoral 
fellow  supported  partially  by  the  NSF,  worked  full  time  on  the  project.  Gregoiie  Winckelmans 
also  worked  full  time  as  a  postdoctoral  fellow  until  August  1993.  Douglas  Shiels,  a  Caltech 
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research  assistant.  Ryan  Mackay,  a  graduate  research  assistant,  spent  the  summer  of  1993  on 
generating  computer  graphics  displays. 
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